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(54) Reduction of image artefacts in 
computed tomography 

(57.) Image reconstruction errors and 
artefacts can occur in computed tomog- 
raphy in which a body section 3 conr 
tained within a positioning zone 25 is 
only scanned completely throughout a . 
smaller scanning zone 2 as a result of 
the incorhpleteness of measuremeht in 
the scanned rejgions 31 and 32 of the 
body section outside the zone 2. 

According to the present invention 
the respective difference value is 
formed between the sum of all absorp- 
tion values measured in a series relat- 



jing to each direction and the maximum. 
! For each direction associated vyith a 
jnon-zero difference value, the measure- 
Imept values are cpmplemented by a 
: priori vaiues bach proportional to the . 
jlength of measurementbeam in the . 
body section but outside the sciannihg . 
zone 2, The boundary of the body sec- 
tion is determined from tangential mea- 
surement beams inside the zone 2 and 
by interpolation therefrom outside that 
zone. A central vector is defined for a 
predetermined set of directions, and the 
boundary is identified as the first radial 
intercept along the vector direction with 
a tangential beam path*. 
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SPECIFICATION 

Method of determining tiie spatial distribution of radiation absorportlon In a body sierctidn 

5 The invention relates to a method of determming the spatial distribution of radiation absorpprtiori in a 5 
planar section of a body from a plurality of measurement series, each of which represents a sequence of 
measurement values, each said value corresponding to the integral of the absorportlon of the body along 
a respective pne of a plurality of parallel measurement bearn paths which pass through a scanning zone, 
the various measurement series being formed from measurement values relating to measurement beam 

10 paths, oriented in corresponding different directions, the planar section of the body to be examined being 10 
completely scanned by corresponding measurement beams when directed in at least one of the direc- 
tions. . M U s 

A method of this kind has been proposed hitherto, for example, from U.K. Patent Specification Number . 
1,283,915. Therein, the absorption distribution is reconstructed from the measurement values by means 

15 of a computer utilizing suitable algorithms. _ 15 

Therein, the scanning zone is a zone which is irradiated throughout and in succession from different 
directions. In the arrangements proposed hitherto for performing the method, the body whose absorption 
distribution is to be determined in one planar section thereof must be positioned insidlB the, normally 
circular, scanning zone. This mfeans that there may be no part of the body section which is not irradiated 

20 by the radiation in any of the measurement directions. OthenA/ise, as tests have shown, senous errors will 20 
tend to occur during the computed reconstruction of the local absorption distribution, even within the 
scanning zone which is itself completely scanned by the measurement beams. These errors may be - 
caused by inaccurate positioning of the body to be examined inside a positioning zone when this zone is 
itself largerthan the scanning zone, or bythefactthatthe body section has an elliptical boundary whose 

25 major axis is longer than the diameter of the scanned zone. ^ , . ^ x • 

Therefore, the present invention has for an object to provide an improved method of the kind referred to in 
which errors caused by the fact that parts of the body section under examination are not irradiated by the 
radiation in at least one direction can be reduced. r u 

According to the invention there is provided a method of determining the spatial distribution of absor- 

30 portion (|t(x,y) of radiation in a planar section of a body from plurality of measurement series, each of 30 
which represents a sequence of measurement values (Q, v, p) each said value corresponding to the integ- 
ral of the absorption of the body along a respective one of a plurality of parallel measurement beam paths 
which pass through a scanning zone in a direction v, each measurement beiam path being a distance p 
from the central beam path of the series,the various measurement series being formedfrom measure- 

35 ment values relating to measurement beam paths oriented in corresponding different directions, the 35 
section of the body to be examined being completely scanned by corresponding measurement beams 
when directed in at least one of the directions, characterized in that the sum of (f(v) of all measurement 

values (Q(v, p)) of a measurement series is formed and that the difference (KM) between each sum and 
the maximum (A) of all sums is determined, and in the case in which at least one of the di^^^ 

40 thus formed departs from zero, a priori wB\ues (Q'*{v,'p)) are formed to complement the measurement. 40 
values, such that for each beam direction V the sum ofthe a pr/or/Values(Q''(v,p)) corresponds to the 

difference (K{p)) determined for the measurement series corresponding to that direction, each a prion 
value (Q'*(v, p)) for each direction being proportional to that length (Mv, p)) of a corresponding beam path 
which extends through the body section in that region which is outside the scanning zone and in s^^^^ 

45 direction, said length corresponding to the distance betweeri the two points of Intersection (x v.Y/i; X2, y2) 45 
of said beam path with the boundary of the body section, said boundary being determined as the en- 
velope of all beam paths which are tangential to the boundary of the body section inside the scanning 

zone, or as an extension of said envelopewhich Is determined by interpolation, the reconstruction of the 
' absorption distribution (|t(x,y)) inside the scanning zone being computed from the measurement values 
5b (Q(v, p)) and the a pr/or/ values (Q^(v, p)). , u- u ^ ^ 

Thus, a pr/or/ values are formed by the method embodying the invention, i.e. values which do not 
follow directly from a measurerhent and which at least approximately represent the absorprtion of the 

body outside the scanning zone, the reconstruction of the absorption distribution inside (and m given cases 

also otitside) the scanning zone is realized on the basis of the measurement values and the a prion values. 
.55 The Idea of replacement of missing measurement values has been proposed (B.E. Oppenheim, "Image 55 

processing for 2-D and 3-D reconstruction from projections", Stanford University, U.S.A., 1975) but no 

suitable method of determining the a pr/or/ values has been described. 
The invention is based on the recognition of the following: . . , 

a; The integral overall measurement or the suni of all the measurement values in one direction; in _ 

60 other words, the sum of all measurement values of a measuring series, forms a value which is substantia - 60 
ly constant and independent of the direction of the beiams in the case of complete scanning, "substantially 
constant and independent" Is to be understood to mean herein that uncertainties of measurement 
caused by the physical nature of the measurement give rise to a small fluctuation of the integral value. 



This fluctuation, which may be the reason why the integral over all measurement values is not exactly 
65 constant when the bodv section is completely scanned; should be.con^idered negligibly small hereinafter.. 



65 
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This assumption is based on the fact that the area integral of the absorption in the body section examined 
is independent of the direction of the pair of rectangular coordinates lying in the body section and along 
which the area integral of the absorption is formed. Because each measurement value represents the line 
integral of the absorption in a given direction, the integral over all measurement values along adjacent 
5 paths each in a given direction i.e. the integral in the direction perpendicular to said given direction is - 5 

equ,at to the area integral of the absorption. Therefore, the sum of all measu^ 
tion, namely of one measurement series, will be a constant and independent of that di^^^ 
. In the case of a part// mcom^/efe scan (during which the body section is dompletely scan 
radiation'in atleast one direction in the section under examination and incompletely scanned in at least one 

10 other direction)/a maximum value will be formed for the integral of the set of measurement values 10 
relating to a direction for which the body islice has been corripletely scanned. In other directions, smaller 
integral values must result. The difference between this maximum Value and the.integral value relating to ; 
given direction will thus correspond to the integral value or the sum of the a pr/or/ values required for that 
direction. ' ^ ■ • : ' 

15 b) Examinations using arrangements of the described kind have shown that the X-ray absorption by 15 
biological tissue varies within only narrow limits in major regions of the human body. Therefore, it may 
be assumed that the absorption inside the body isconstantfor a given direction along all beams situated, 
outside the scanning zone, 
c) The outer boundary of a section of a human torso or head may be represented as a.clpsed and 

20 shiooth curve. For regions of the body which are situated inside or at the edge of the scanning zone, the 20 
position of the outer boundary can be derived directly from the measurement values as an envelope 
formed from all measurement beam paths which are tangential to the boundary of the body section. For 
those parts of the body section which remain outside the scanning zone, the position of the boundary 
curve can be approximated by interpolation. The length of the beam paths passing through the body 

25 section which are assumed to be situated outside the scanning zone is thus approximately defined. / 25 
The a ;pr/or/ values can be calculated on the basis thereof. 

An embodiment of the invention will now be described by way of example, with reference to the accom- 
panying drawings, of which:— 

Figure 7 shows the geometrical arrangement of a section of a body to be examined with respect to the 
30 scanning zone and the positioning zone, - 30 

F/fiTL/res 2a and b illustrate the variation of the measurement values of a measurement series, ' 

Figures 2c and d respectively show the sum value of the measurement series as a function of the 
irradiation direction for a completely scanned body section and fpr a partly incompletely scanned body 
section, 

35 Figure ^3 is a diagram which shows an example of the determination of a central vector, 35 
Figure 4 shows the variation of the central vectors as a function of the polar angle <p. 
Figure 5 sho\A/s diagrammatically an arrangement for performing a method embodying the invention. 
Figure 6 illustrates the variation in the measurement values and in the a pr/or/ values, and 
Figure 7 is a diagram which further illustrates the geometrical arrangement of an incompletely scanned 

40 body section: . y' 40 

: Figure 1 illustrates a circular scanning zone 2, having a centre 1, which is generally also the centre of 
rotation of an assembly formed by a radiator and a detector array. The radiator, the detector array, the 
rnotors and the drives and all other mechanical parts are not shown since they may be of similar construc- 
tion to that employed in apparatus for computed tomography hitherto proposed (for example, as in the 

45 apparatus described in U:IC Patent Specification No. 1,283,91 5). Regions 31 and 32 of a body section 3 to . 45. 
be exarnined, are situated outside the scanning zone 2, but inside a positioning zone 25 which is assumed . 
to be circular and concentric with the scanning zone 2. The outer boundary of the body shown in Fig, 1 
corresponds to a section through the upper thorax and the shoulder bones. The measurement beam 
paths 4 and 5 (which are the beams or strip paths along each of which the absorption is determined and 

50 represented by a measurement value during a measurement), belonging to two different series, are 50 
shown by way of example, said beams irradiating the scanning zone 2 along different angular directions 
but always scanning the whole of the scanning zone 2. The central beams passing through the centre 1 of 
the scanning zone 2 are denoted by the references 41 and 51 , while the other measurement beams which 
are tangential to the scanning zone 2 are denoted by the references 42, 43 and 52, 53, respectively. It can 

55 be seen that the measurement series denoted by the reference 5 completely covers the body 3 and that 55 

the measurement Series 4 does not completely cover the body 3, the shaded areas 31 and 32 not being 
. covered. ' . . , 

The Rgures 2a and 2b are graphs 54 and 44, respectively, of th measurement values Q(v, p). measured in 
the measurement series 5 and 4. Q(p, p) serves to express that each measurement value is a function of 

60 the direction measured from the vertical in the figures, in which the measurement beam irradiates the 60 
body during the determination of the measurement value, this direction being the same for al l measure- 
mient values of a measurement series, and also a function of the position p of a measurement beam path 
with respect to a straight line extending perpendicularly to the measurement beam path direction, and is 
the distance between the relevant measurement beam and a corresponding beam through the scan cen- 

65 tre, the said distance p relating to measurement beam paths which are situated to one and to the other 65 
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side of the central beam path, having opposite signs. It will be apparent from Figures 2a and 2b that the 
curve 54 varies relatively gradually, while the curve 44 goes abruptly to the value zero atthe edges 45 and 
46 where it represents the measurement values along the outer measurement beam paths 42 and 43. 
Each of the measurement values represented by one of the curves 44 and 45 can be calculated in known 
5 linanher by means of the formula: 

= J^y" (^ty)^s . ; (1) . 



10 



15 



20 



Xt/':) = J Q(l^., P)clp = JjT /u(x,y)dp ds (2) 



10. 



in which Q{v,pHs the measurement value determined along the measurement beam path having a direc- 
tion V and a distance p from the scan centre, and ji(x,y) is the absorption of the body slice at the point (x,yjf 
of a fixed rectangular system of coordinates. The quantity sis a coordinate extending along.the direction 
of the rheasurement beam path. . 15 

When the integral f (v) over all measurement values of a measurement series is formed, the following is 
obtained: . ' 



20 



in which p and sform a coordinate system rotated in synchronism with v and having a fixed coordinate 
25 origin at the scan centre 1 . Because this integral f (v) corresponds to the area integral of the absorption of 25 
the body over the plane of examination, which Is independent of the direction v, as has already been - 
stated, this independence must also be applicable to f(v), so that for a complete scan of the body section 
the graph 6 of the function f(i/) in dependence on v results, LeA(v) is independent of v and is of a constant 
value = A. This relation is applicable whenever the measurement value Q (v, p) is proportional to the line 
30 integral of the absorption in the direction v along the measurement beam path having the positional 30 
distance p. This property, assumed to be a given factor hereinafter, can be obtained in known manner, for 
example, by forming the logarithm of the intensity ratio of the unattenuated measurement beam to the 
measurement beam attenuated by the body. Furthermore, the integral f(v) is not formed by a real mtegra- 
tioain practice, but rather by summing all the measurement values in a measurement series. However, 
35 this does not affect the validity of the considerations underlying the present embodiment. .: 

For a measurement series which does not completely cover the body section 3, the value of f(v) must be 

less than this constant value: 



40 



^<Z/> - / Q(^,p)dp < A (3) 



50 



A s= Max 



40 



Therein, R is the radius of the scanning zone and A is the previously mentioned constant value. During an 
45 examination involving a partly incomplete scan of the body section (i.e. during an examination in which 45 
the body is completely scanned during some of the measurement series, whilst it is not completely scan- 
ned during other said series), the constant value A can be obtained as a maximum of the various values 
f{v): , 



50 



55 Fig. 2d shows, byway of example, a graph 65 of the integral valuef(v) forthe case (shown in Fig. 1) 55. 

concerning a partly incomplete scan, the references 47 and 57 denoting the direction of the measurement 

beams of the measurement series 4 and 5, respectively. 
The equations (3) and {3a) are used in the embodiment of the invention to establish whether the case of 

a partly incomplete scan.has occurred and for calculating the integral value K(v) of a pr/or/ values Q"" (v, P) 
60 necessary to complement the measurement values. This. is because: . 6p 



65 



. - . 65 



GB 2 009 560 A 



If K (v)is greater than zero, this is proof pf an incomplete scan only if no errors have occurred in the 
various rneasurernent series. On the other haridjn the case where it has been ascertained th 
has been completely scanned, an error (caused, for example, by movements of the body or fluctuations of 
the detector sensitivity) may be assumed when K (v) is greater than zero. 
5 When it can be assumed^ as has already been mentioned, that the absorption at or adjacent the bound* . 5 
. ary of the body section, for the most part departs only slightly from a mean value, the equation (4) can be 
. used for determining the a pr/ori values in the regions 31 and 32 (Rg. 1) outside the scanning zone 2, 

provided hat the area of the body section outside the zone Is known. In that case, the following must be 
. applicable for each a pr/on value Q^(v, p): 

10 \ * 10 

. Q*(2/ ,p) = c(i/^) . i.{2/,p) «5) 

15 15 

wherein C(v) is a proportionality factor and L(v, p) is the length, yet to be determined, of an (assumed) 
beam path which irradiates the body in that region lying outside the scanning zone but inside the body. - 
The following is then obtained in cbnibination with the equation (4(; 

20 r ^ ^ 1 26 

25 Therefore, when the lengths L(v, p) have been determined, which will be described hereinafter, the 25 
factors C (v), being constant for a given direction v, and hence the a pr/or/ values Ol^{v, p), can bd calcu- 
lated by means of the relations (4), (6) and (7). Forthe determination of L(y, p), the embodiment utilizes the 
fact that the body section generally has a closed, smooth and substantially convex boundary. This bound- 
ary can be described by a series of central vectors r (9), where: 

30 : , ; . 30 



f)^ y^i^) (8) 



35 : ■ 35 

where x(if>),y(<p) is the point of intersection of a central vector, originating from the centre 1 of the scan- 
ning zone 2 in the direction <p, with the boundary of the body section. When the vector lengths lr(9)| are 
plotted as a function of 9, a substantially gradual curve is obtained. 
For determining the boundary of the body section, firstly those beam paths p3 and p4 (see the beams 1 1 4 
40. ... 1 17 in Fig. 3) must be determined which are tangefntial to the boundary of the body section under 40 
examination. Forthese'measurement beam paths, the following is applicable in the interval between p^ ^ 
and Pa: ■ 

(i(^, p)>Oor V . / 

Q('&, p) = 0 outside the interval [pa/pj (9) 
45 Fig. 3 shows the points of intersection of a radial line 1 30, extending in the direction 9, with tangential .45 
bearii paths 114, 115 and 1 16 thus established,. said points being denoted by the references 120, 121 and 
: 122. It can be seen that the point of intersection 122, being situated at the smallest distance along the 
radial line. 130 from the centre 1, substantially coincides with the edge of the body-This relation can be 
generalized: 

50 * \ . - . 50 



7) 



55 ■ - . • . . 55 

in which x((p, ^) and y(cp, ^) is the point of intersection of the radial line extending in the direction 9 and a 
beam which extends in the direction ^ and is tangential the boundary of the body section; for x(<p, -&) and 
y(<P,^): 

x(<p, d) = P3,4('&) • 1/(sin tan 9 + cos ^) ,^ ^ , 
60 y(<p, •&) = p3/4('&) l/(sin + cos ^ - tan 9) 60 

wherein the angle in this expression is the.angle between a straight line which extends horizontally to the 

right (in Fig. 3) from the centre i and the perpendicular from the centre to the tangentially extending 

measurement beam path. 

In order to prevent successive errors during the determination of the vectors r(9) in the case of a beam 
65 path direction for which the positioned distance of pa = R and/or of p4 = — and K(^) > 0, the relation (9) . gg 
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must be completed by , 

or (9a) 

P4 = -^Ro 

5 Where Rois the radius of the positioning zone 25 (see Fig. 1 ) and the sign of Ro depends on whether the 5 

position of p3 or P4 is to one side or to the other with respect to the central beam. For the conf iguration 
. shown In Fig. 3, the variation of the central vectors r(<p) as a function of <p as shown in Figu4 in thus 
obtained. The reference of 104 denotes that portion of vectors rt^p) which correspond to the boundary region 
of the body section lying outside the scanning zone 2 to the left in Fig. 3. In accordance with the equations 
10 (9a), (10) and (1 i), the maximum vector length Is limited to approximately Ro- The reference 105 denotes 10 
that portion of the vectors r(((>) which correspond to the boundary region of the body section lying outside 
' the scanning zone 2 to the right in Fig. 3. Because the point of intersection 1081 (Fig. 3) is in this case 
situated inside the boundary 25 of the positioning zone in the device, the vector lengths are not limited. 
The references 1041 and 1051 in Fig. 4 represent the variation of the true vector lengths Ir(cp)| in accord- 
15 ance with the relation (8). \. _j - *- \f / V 

In order to reduce the difference between the two curves, and hence the errors m the determination or r(<pj, 
the vector length determined In accordance with the equation (10) is replaced by a corrected value (<p), 

where: _ ; — . , . . - - 

?* (cpr= r(9i) ' ai + r(<p-i)a2 + r(<p2) - as + r(cp -2) a« + ... (10a) . . 
20 The interpolation iri accordance with the equation (iOa) is performed whenever |r(<p)|* determined in . 20 
accordance with the equation (10), has a value which is greater than the radius R of the scanning zone. In 
this respect, use is made of values r((pi) . . . r((p2) which are smaller than or at the most equal to R, their 
angular.positions 91 . . . <p? being nearest to the position 9-. The terms ai, 32, as, 34 . . . are interpolation 
coefficients which are fixed in accordance with known interpolation methods, and the angles 9satisfy the 
25 relationship: 9 2 < 9 -i < 9 ,< 9 2. The boundary of the body 3^^^ 25 
vectors rCcjp) or r^(9). *. 

Using r(9) and r^(9), each one of a set of lengths U-^, p) can be calculated as the distance between two 
corresponding points of intersection (xi, yO and (X2;y2) of a beam path which extends outside the^canning 
zone 2 in a direction ^ and having the positional distance p > R, with the boundary of the body section 3. In 
30 this respect: 



30 



35 



35 



(12) 



where Xi = |r(9')l*cos9',X2 = |r{9")l- cos 9" 

yt = fr(9')|.sin 9^y2 = |r(9")l. sin 9" x - 

The polar angles 9' and 9" of the two central vectors, terminating in thesaid points of intersection, are 
40 determined jn accordance with the following condition: ^ 

|p/cos(f^-9')"lr(cp')J< E ^^gj 



or 



|p/cos(^-9")-|r(9") |< 2. . : , , C^R/mn 

Therein,gis a small value in comparison with the radius R of the scanning zone, for example, R/100. 



45 



45 The arigle 9 is then varied until the relation (13) is satisfied 

Hereinafter, the calculation process. will be described for performing a method.embodying the >nyen- 
tioh by means of a computer which is diagrammatically shown in Fig. 5. The measurement values Q(^, p) 
arrive from a memory (not shown), via a data line Vm anarjthmetic unit 2 in which the integral or the sum 
of all measurement values f(^) is formed in accordance with the equation (3) for each measurement series, 
50 /-e. for each The Vjalues f(^) thus calculated are stored in a memory 3. An arithmetic unit 4 determines, in 50 
accordance with the equation {3a), the maximum value A of the values stored in the memory 3. An arithmetic 
unit 5 determines, on the basis of the values f(^) stored in the memory 3 and the maximumm value A 
calculated in the arithmetic unit 4, the differences K (^) in accordance with the equation (4) and stores these 
values in a memory 6. From the mehnory 6 a control logic line 7 extends to a switch 8 which enables a signal 
55 path 9 for the measurement values Q('&,p) when all the differences K(^) are zero or less than a . ; 55 

predetermined low threshold value. The measurementyalues Q(^, p) can then be applied, without further 
modification, via a data line 10, to a known arithmetic circuit (not shown) for determining the absorption. At - 
the same time, all further calculations in arithnrietic units 12 to 24 are inhibited via a control logic line 1 1- 
However, in the case of a partly incomplete scan, so that at least one value K{^) > 0 exists, the switch 81 is 
60 closed via the control logic line 1 1 , so that the measurement values Q (^, p) can be ajDplied to the arithmetic 60 
unit 1 2 w'a a data line 91 . This unit then determines, for all directions ^, the positions of the 
(measurement) beams which are tangential to the body by means of the equations (9) and (9a) and stores 
this data in a memory 13. The points of intersection x(9, d), y (9, ^) for all beam directions ^ and, prefer- 
ably, 360 equidistant directions 9 of the central vector, are determined in an arithmetic unit 14 in accord- 
65 ancewiththerelatibn (ll). An arithnnetic unit 15 then determines, by means of the relation (10), the 65 
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lengths |r(<p)| and stores this data in a memory 16, The interpolation in accordance with the relation (10a) 
is performed in an arithmetic unit 17. If the value K (^) becomes greater than zero for a given beam direction 
^ (which impiies that the body was lipt confipletely scanned in the direction d), the completion of 
the calculation process in an arithmetic unit 19 is enabled via the control logic line 18, the unit 19 calculat- 
5 ing, using the equation (13), for all beam positions |p|^Ro, the central vectors r(q>') and r(M>") whose end 5 
points (xv y,) and (xz, Vz), respectively, are situated onthe beam path relating to the position p. If such end 
points exist, Le. if the beam path relating to the position p intersects the boundary of the body section, the 
associated length U-a, p) is calculated in an arithmeticunit20in accordancewith the equation (12). 
For beams which are situated outside the scanning zone 2 and {Rg. 1 ) and outside the boundary of the 
10 body section, and also for all beam directions ^for which the calculation is not enabled in the arithmetic 10 
unit 19 (because the body \s completely scanned in the direction ^, so that K(^) = 0), the arithmetic unit 20 
produces L(^, p) =^ 0/ AH values L(-&, p) are stored in a memory 21 . The values K(^) made aya^^^^ 
data line 23 are used in an arithmetic unit 22, in conjunction with the values L(-^. p) stored in the memory 
21 , for determining the factor G{^) in accordance with the relation (7). On the basis thereof^ the arithmetic 
-J 5 unit 24 determines the a pr/on values Q'' p) to be complemented and applies theseyalues via 10 to the 15 
arithmetic unit (not shown) for determining the absorption distribution. . 

Rnally, it is to be noted that in the case of incomplete scanning in all beam directions Le. in the case 
where a measurement value departing from zero is determined for all beam directions at least one 
portion of the boundary of the scanning zone (QC-^ R), or Q(^, -R) > 0), the arithmetic unit 12 provides an 
26 indication, via a control logic line 25, on a data display (not shown) that a satisfactory determination of the 20 
absorption is not possible in these circumstances. The arrangement shown in Fig. 5 can be realized by 
'means of known digital computer components. 

The method embodying the invention can bie further refined. For example. Fig. 6a shows a measure- 
ment value profile complemented by a pnon values Q''(^, p) 451 and 461. As is clearly shown, at the 

25 boundaries 45 and 46 between the curve 441 which represents the measurernent values Q(^, p) and the 25 
curves 451 or461 which representthe a pr/onvalues d'^i-^, p) an irregularity is present which rtiay be 

' . caused, for example, because the circuiar scanning zone 2 has intersected a part of the body having a 
higher absorption. In order to diminish such irregularities and to endeavour to ensure that a curve of the 
form indicated by 452 (Fig: 6b) is pi^ovided, the values Q*'(^, p) are corrected by interpolation in an im- 

30 proved arrangement to provide a modified a pr/on value given by:~ 30 
QXX p) ^ ^ Q R) . bp + Q'^ p) . C2 {^) (14) 
• where bp is a weighting factor for which:* 

bp = 1 for p = R and bp = 0 for p > R and for Q"* p) = 0. 

The f actor C2(^) is chosen so that the area integral under the curve 451, representing the values □*'(•&, p), is . 

35 equal to the area integral under the curve 452 which represents the corrected value Q''"'!^, p). As a result, it 35 
can be arranged that the integral over all values associated with a given beam direction ^ is constant ( = 
the maximum value A). The value 

AQ(^,R) = Q('&,R)-Q^(^,R)C2(^^) (15) . o 

is defined as the difference between the measurement value Q (^, R) at the edge of the scannmjg zone 2 

40 and the modified a pr/or/ value Q'*(^, R) along the measurement beam path for which p = R and is added, 4(0 , 
weighted by bp, to the a pr/ori values in order to ensure a smooth transition at the point where p R (45 in 
Fig.6). u 

Afurther improvement can be achieved when the vectors ri^p) determined in accordance with the equa- 
tion (10) are taken into account, by contrast with the conditions described with reference to the equation 

45 (10a), also for calculating the edge of the body in accordance with the equations (12) and (13) when |r (tp)I 45 
>R, subject to conditions which will be explained hereinafter. As a result, the number of vectors which are 
to be replaced by interpolated values r' dp) is reduced, so that the edge of the body can be more accurately 

determined. « u 1 • 

For an explanation of the required conditions, the beam positions calculated in accordance with relation 

50 (9) are shown in Fig. 7 for various beam directions Therein, 531 and 532 denote groups of beams which 50 
extend through the scanning zone 2, and the reference 431 denotes a group of beams which extends 
outside the scanning zone 2 and for which p is made equal to ±Ro in accordance with the equation (9a), 
the length of the vectors 61 and 62 corresponding exactly to the radius R of the scanning zone 2. It can be 
seen that there are also vectors whose length exceeds the radius of the scanning zone 2, for example, the 

55 vectors 63 and 64 which correctly describe the boundary of the body section. However, not all vectors - 55 

between the vectors .63 and 64 are suitable for describing the boundary of the body section 3, so that they 
must'be replaced by interpolation in accordance with the equation (1 Oa). In order to deternriine in which 
the angular region of the central vectors must be replaced by interpolation^ use can be made of the fact 
that a gradually monotonously decreasing or increasing reliatiqnship exists between the angle 9 of the 

60 central vector arid the angled of the measurement beam path on which the central vector terminates for a 60 
body section whose boundary is convex in the plane of examination. Therefore, when the determined 
' central vector equals the radius Ro of the positioning zone 25, or when the central vector terminates each, 
time on the same measurement beam path (p = R) extending from the outer edge of the scanning zone 2, 
the value of R(<p) must be determined by ihteipolation in accordance With the equation (10a). 

65 The method embodying the invention can also be used for scanning systems in which a group of 65 
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diverging beams simultaneously completely scans the body section, because groups of beams of this 
kind can be re-arranged to form groupsof parallel beams by known methods. Scanning systems of this 
kind are sometimes referred to as third-generation scanners; the scanning zone is then completely irradi- 
ated by a fan-shaped radiation beam which originates from the radiator. Similarly, the method embody- 
5 ing the invention can also be used for so-called second-generation scanners. It is also to be noted that the 
method embodying the invention of X-radiation, but also for the determination of the distribution of the 
emission of gamma-radiation utiliiirig a spri^lled radionuclide scanner. 

^'l^"^method of determining the spatial distribution of absorption (n(x, y)) of radiation jn a planar 

10 section of a body from a plurality of measurement series, each of which represents a sequence of mea- 
surement values (Qlv, p)) each said value corresponding to the integral of the absorption of ths body 
along a respective one of a plurality of parallel measurement beam paths which pass through a scanning 
zone in a direction v, each measurement beam path being a distance pfrom the central beam path ottne 
series, the various measurementseries being formed from measurernent values relating to measurement 

15 beam paths oriented in corresponding different directions, the section of the body to be examined being 
completely scanned by corresponding measurement beams when directed in at least one of the direc- 
tions, characterized in that the sum (fu)) of all measurement values iQ{v, p)) of a measurement series is 
formed and that the difference |K( v)) between each sum and the maximum (A) of all sums is determined, 
and in the case in which at least one of the differences (K(y)) thus formed departs from zero, g pr/oc/ values 

20 (Q'-Cv p)) are formed to complement the measurement values, such that for each beam direction v the sum 
of the a pr/or/ values (Q^v, p)) corresponds to the difference (K(p)) determined for the measurement series 
corresponding to that direction each a pr/or/value (Q'<(i', p)) for each direction being proportional to that 
length (Uv, p)) of a corresponding beam path which extends through the body section in that;regipn 
which is outside the scanning zone and in said direction, said length corresponding to the^distance be- 

25 tween the two points of intersection (x„ y,; x.,y.) of said beam path with the boundary of the body section, 

said boundary being determined as the envelope of all beam paths which are tangential tothe boundary 
. ofthebodysectioninsidethescanningzone,orasanextensiohpfsaidenvelopewhichisdeterminedby 
' interpolation; the reconstruction of the absorption distribution (,i.<x,y)) inside the scanning zone being 

computed from the measurement values (Q(v, p)) and the a pr/or/ values (Q''(v,p». 
30 2 A method as claimed in Claim 1, characterized in thatforthe determination of the measurement 

beam paths which are tangential to the boundary of the body section, those measurement values are 
determined which, when viewed from the centre of the scanning zone, first produce approximately the 
value zero, whilst all parallel beam paths which are situated further outwards from said centre also pro- 
duce the value zero, the identities of the position {p.. p.) or the respective distances p of the measurement 
35 beam paths associated with these values, relative to. the centre of the scanning zone, being stored. 
3 A method as claimed in Claim 1, characterized in thatwhen a measurement value which departs 

from zero is determined for a measurement direction (v) at the outer edge of the scanning zone, the radius 
(R.) of a positioning zone largerthan the scanning zone and which encircles the body section, is set as a 
position (p3 or p4 of the beam path tangential to the boundary of the body section 40 
40 4 A method as claimed in Claim 1, characterized in that the envelope is formed as the locus of the 
vector termini! of a progressive series of central vectors M^))) each pfwhicKconnects the centre of the _ 
scanning zone to the respective points of Intersection of a radial line extending along the direction of said 
central vector with that beam path tangential to the boundary of the body section vyhich IS situated at the 

shortest distance along said radial line from the centre of the scanning zone. - 45 
45 5 A method as Claimed in Claim 4, characterized in that the central vectors are formed for a finite 
number of angular positions (9) which are uniformly distributed over one complete rotation 

6 A method as claimed in the Claims 1 to 4, characterized in that those central vectors (r (<p)) yvhose 
lengths are equal to or greater than the radius (R) of the scanning zone, are replaced by corrected central 
vectors (r" (<p)) which are determined by interpolation from the nearest adjacentcentral vectors (r (<p..) .... ^ 

50 rivt)) whose lendths are shorter than the radius (R) of the scanning system. ' . ' , . / 

• 7 AmethodasclaimedinClaim1,characterizedinthatineachofthemeasurementdirecfions(i.)for: 

which K(v) is not zero, further beam paths ara defined which are situated outside the scanning zone but 
inside a positionirig zone which encircles said scanning zone, and corresponding pairs of central vectors (r 
(<p'), r(<p")) are determined whose vector terminii are situated approximately on these defined further beam 

55 paths, the distance between the two vector terminii being used as a measure for the a prion value (Q )w. 9 
p)) of those beams to be complemented. j. . , ,»u • ♦ 

8 A method as claimed in Claim 7, characterized in that in each measurement direction M the integral 
is formed overthe distances along all those beam paths which are situated outsidethe scanning zone but 
inside the positioning zone, a proportionally factor (C(v)) being formed from the ratio of the difference 

60 (K(v)) to this integral, said proportionally factoi'fornfiing the a priori value (Q>') after multiplication by the 
distance (L(v, p)) between the two Vector terminii (r (ip'K r(<p"). - , 

9. A method of determining the spatial distribution of absorption of radiation m a planar section of a 

body substantially as herein described with reference to the accompanying drawings. " . 
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